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Abstract 

We present a new Fortran code to calculate the scalar one-loop four-point integral 
with complex internal masses, based on the method of 't Hooft and Veltman. The 
code is applicable when the external momenta fulfill a certain physical condition. 
In particular it holds if one of the external momenta or a sum of them is timelike 
or lightlike and therefore covers all physical processes at colliders. All the special 
cases related to massless external particles are treated separately. Some technical 
issues related to numerical evaluation and Landau singularities are discussed. 



1 Introduction 



In order to calculate the radiative corrections to reactions with unstable internal particles 
(like the top quark or W/Z gauge bosons) which can be on-shell, one has to resum the 
propagator of these particles. This can lead to prescriptions where the squared mass 
of unstable particles are regarded as complex parameters with a finite imaginary part 
proportional to the width. It means that we have to calculate Feynman integrals with 
complex masses. 

From the mathematical viewpoint, the introduction of a width helps to protect the 
cross section from several singularities, called Landau singularities [H El El IH El |6] . Let 
us explain this in a simple way. In order to define a Feynman integral in the case of 
real masses, one needs the ie-prescription to deform the integration contour from various 
poles. The final result is obtained by taking the limit e 0. We observe that the position 
of the poles, in general, depend on the values of internal masses and external momenta 
(see Appendix |A] for more details). It is therefore possible that these poles can pinch the 
integration contour in the limit e ^ 0, leading to a Landau singularity. By making the 
masses complex, the poles are moved away from the real axis and hence all the various 
Landau singularities no longer exist. 

The method to calculate scalar one-loop integrals with real/complex masses was ex- 
plained by 't Hooft and Veltman [7] . Following this, van Oldenborgh has written a Fortran 
package FF to calculate one-loop integrals [HI Ej . Based on this package, a more convenient 
library called LoopTools which can calculate scalar and tensor one- loop integrals up to 
five-point has been developed by Hahn and colleagues [IHllII]- LoopTools/FF can handle 
one- loop integrals with complex masses up to three-point. The only missing piece in the 
scalar integral sector is the one-loop four-point function with complex masses. A Fortran 
77 code to calculate this function is presented in this paper. The method is based on 
the one given in [7] where a general result is written down, see Eq. (6.26) of [7j. This 
result is applicable as long as the external momenta fulfill a certain physical condition. In 
particular it holds if one of the external momenta or a sum of them is timelik^. If there is 
at least one lightlike external momentum then similar results can be derived by using the 
same method. If the condition is not met then the calculation becomes mathematically 
quite tricky. In principle, the result in this case should be an analytical continuation of 
Eq. (6.26) of [7] with some extra logarithms. However, it is not easy to find these terms in 
practice. Fortunately, this special case is never encountered in collider scatterings where 
there always exists at least two timelike/lightlike external momenta. 

The outline of this paper is as follows. In the next section, we give all the relevant 
formulae which have been implemented in our code, named DOC. All the important facts 
related to the code are given in Section El The appendices include a discussion on three- 
point and four-point Landau singularities and the calculation of a basic integral. 

*Iii our metric timelike p implies positive p^. 
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2 The method 



The notation used in this paper follows closely the one of Ref. [7j. The definition of 
logarithm, Spence function and 77-function are the same. The signature of our space-time 
metric is (+, — , — , — ) which is different from the one used by 't Hooft and Veltman, 
however. 




P: 



Figure 1: The box diagram. 



With the momenta assigned as in Fig. [H the scalar one-loop four-point function in 
4-dimension space time is defined by 



Do 



1 



2 — m'j + ie] 



where the momenta fcj are related to the real external momenta as 

fcl = 0, k2=Pl, k3=Pi+P2, k4=pi+p2+P3 

and the internal masses can have a non-positive imaginary part 

= mf — irriiTi, mjFj > 0, z = 1, 2, 3, 4 



(2) 



(3) 



with Ti are the width of the internal particles. The Feynman parameter representation 
reads 



Do 







°° 5(y]^_ Xi - 1) 

dxidx2dx3dx4 ^~^2^ ' 



A = - ^ XiXjQij - ie, Qij = m} + m^^ - {ki - kj 



(4) 



2 



Qij is called the Landau matrix. By changing the integration variables as t = Ylt=i 



X. 
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X 



Si=i Xi, y = Xi + X2, z = Xi and integrating out the Dirac delta function we obtain 



pI px py 

Dq = dx dy dz 
Jo Jo Jo 



{ax"^ + by'^ + gz"^ + cxy + hxz + jyz + dx + ey + kz + f — ieY 



(5) 



with 

a = pli b = pl, g = pl, c = 2p2-P3, h = 2pi.ps, j = 2pi.p2, 

d = ml-ml-pl, e = ml - ml - pi- 2p2-P3, k = ml-ml+pl + 2pi.pi, (6) 
/ = ml. 

It is important to notice that a, b, g, c, h and j are real while d, e, k and / are complex. 
The denominator always has a negative imaginary part as indicated by the " — ie" 

Im (dx + ey + kz + f — ie) = Im [zm\ + {y — z)m\ + [x — y)m\ + (1 — x)m\ — it] < 0. (7) 

The Landau determinant as a function of new variables reads 

det(g) = k'^c^ + e^h^ + d^f - 2ehkc - 2djkc - 2dehj 

+ Acdeg + ^cfhj - Afgc^ - Aaff - Aae^g + Aaejk (8) 
- 4bgd^ + Abhkd - Abfh^ - 4bak'^ + IGbafg. 

The scalar one-loop four-point integral is a function of 10 variables: 

Do = Do{ai, 02, CIS, ^4, as, ae] bi, 62, &3, ^4) (9) 

where the set of variables are given in Table [1] which shows all the symmetries of this 
functiorQ. Using those symmetries, we see that the notation of external momentum is 
extended to include also (pi + P2) and {p2 +^3)- This definition of external momentum 
will be used hereafter. 



2.1 At least one light like momentum 

If there exists at least one lightlike momentum which we can always choose to be pi then 
the denominator in Eq. ([5]) becomes linear in z. The last integral can be easily performed 
to get 

Do^ = '^{a,b,c,hj;d,ej,k) - T{a,b + j,c + h,hj;d,e + kj,k), (10) 

where the calculation of the T function is given in Subsection 12. 3[ The result contains 72 
Spence functions since each T function contains 36 Spence functions. 

In the case where Dq has two non-adjacent lightlike momenta, say pf = pf = 0, the 
result can be written in terms of two T functions as above. However, one should rather 

tEq. (HI shows that Dq has 4! = 24 symmetries. 
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Table 1: Symmetries of Dq 
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combine them into one single function as (in this way, the result contains a lesser number 
of Spence functions) 

D^o^^ = dx dy 
Jo Jo 

y 

{Gy"^ + Hxy + Dx + ly + F - ie){Ay'^ + Cxy + Dx + Ey + F - ie) ^ ' 

where 

A^h + j, C ^c + h, D ^d, E^e + k, 

F^f, G^b, H^c, 7 = e. (12) 

We notice that each denominator is a linear function in x hence one can perform the 
x-integral by using: 

dx dy = dy dx. (13) 
Jo Jo Jy 



The result is written in terms of 32 Spence functions [3l [6] 



SV -TU \ U + V Uy + V 

2 4 



i=l j=l 



where 



S = Hy + D, T = Gy'^ + Iy + F -ie, 

U = Cy + D, V = Ay"^ + Ey + F -ie; (15) 
and 1/1,2 are the two roots of the equation 

{HA - CG)y'^ + {AD + HE-DG- CI)y + DE - DI + {H - C){F - it) = Q (16) 

oi@ 

_-{AD^HE-DG-CI)T\/^^{^ 
y^'' - 2{HA - GG) ^^^^ 

where the indices 1, 2 correspond to — and + signs respectively and the discriminant, 

det(Q) = {AD + HE-DG- Gif - A{HA - GG){DE -DI + HF - GF), (18) 

is nothing but the determinant of the Landau matrix Q defined in Eq. ([8]). The various 
coefficients of the argument of the logarithms are given in Table [2l The integral in Eq. 
will be considered in Appendix [Bl 



Table 2: The coefficients of the Dq -function 

j Cj^ 

1 A E + G D + F-ie 

2 G I + H D + F-ie 

3 G + H D + I F-ie 

4 A + G E + D F-ie 



2.2 General case 

As proven in (see Eq. (6.18)), Dq can be written as 

Do = —T{a + h + c,g,i + h,c + 2h + ha + ja, j + 2ag; d + e,k, f,e + ka) 

+ (1 — a)T{a, b + g + j,c + h,c + ha, {j + 2ag){l — a)]d,e + k, f,e + ka) 
+ aT{a,b,c,c + ha, —{j + 2ag)a; d,e, f,e + ka) (19) 

■fin the case of real masses, the "je"-part of each root should be written separately as shown in 
Eq. (TOll . 
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where a is a root of the equation 



ga^ + ja + b = 0. (20) 

The method is restricted to real a, i.e. 

f-Abg = A[{pi ■ P2Y - pIpI] 

= \[{Pi+P2f,pl,pl]>0, (21) 
X{x,y,z) = + y'^ + -2xy -2xz -2yz. (22) 

This is satisfied if pi(or p2 or pi + P2) is timehke. If the above condition is not fulfilled 
then one has to keep rotating the external momenta by using the symmetric properties 
given in Table [T] until it is met. The right hand side (rhs) of Eq. f|T9|) contains 36 x 3 = 108 
Spence functions. 

2.3 The r function 

We now consider the following T function 



r(A, 5, C, G, H] D, E, F,J)= [ dx [ 

Jo Jo 



dy 



1 

^ {Gx + Hy + J){Ax^ + %2 + cxy + Dx + Ey + F - le)' ^^^^ 

where A, B, C, G and H are real and D, E, F and J can be complex with the restriction 
that the imaginary part of the second denominator has a definite sign, say negative as 
indicated by the —ie. In order to make the imaginary part of the first denominator 
also negative, we multiply both the numerator and denominator of the integrand with 
—sj = — sign(J) and get 

'^^~^'^Io L '^^ {Gx + Hy + J- t€') {Ax^ + By^ + Cxy + Dx + Ey + F - te) ^^^^ 
where G = —sjG, H = —s jH, J = —sjJ and e' is infinitesimally positive. 

2.3.1 General case: A^O and B ^0 

Following the method explained in we obtain (similar to Eq. (6.26) of [7]) 

with the following cascade of notation. 

-L- v/det(g) -L + v/det(g) 
yi = ^ , Z/2 = ^ , det(Q) = L - ANK (26) 
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are the roots and the discriminant of the equation 



Ky^ + Ly + N 
L 
N 



0, K = B{G + pH) - H{C + 2PB), 

{G + (3H)E - H{D + (3E) - J{C + 2/35), 

(G + f3H){F - ze) -{D + (3E)J 



with P is one root of the equation 

Bf3^ + Cp + A = 0. 



(27) 



(2^ 



It is easy to check that the discriminant is free of (3 and is indeed the determinant of the 
Landau matrix Q if our calculation is set in the context of Eq. (ITOi) or Eq. (|T9l) . And 



Ri 
Pi 



G + pH 



Byf + Ey, + F-te {G + 2pB)y, + D + PE' 

3 .1 



dy 



c 



■J 







i-iy+^Cjy-bj-y^ 



1 



dy- 



(29) 



where A/j is the beta triangle [0,1 — P, —P] as shown in Fig. [2l The second term in 
Eq. fl25l) is nothing but the residue contribution. All the other coefficients are listed in 



-3 


Im 


1- 


\ 


y. 





Re 



Figure 2: The beta triangle [0, 1 — P, —P] in the complex plane. 



Table El 



Table 3: The coefficients of the T function 



j 


9 


6, 


A, 




G, 


1 


1 


P 





H 


G + J-ie' 


2 


-1+P 








G + H 


J -It' 


3 


-P 








G 


J -ie' 


4 


-1 


P 


B 


G + E 


A + D + F-ie 


5 


l-P 





A + B + G 


E + D 


F -ie 


6 


P 





A 


D 


F -ie 
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We have an important remark here. Eq. (125!) is very different from Eq. (6.26) of ^ 
in the respect of numerical evaluation, 't Hooft and Veltman suggested a way to proceed 
with the latter by splitting the residue into two parts: one part has the same form as the 
first term in Eq. fl25l) but with the variable y in the argument of the logarithms replaced by 
the pole of the corresponding previous factor, the other part is related to the r^-function. 
Eq. (6.26) of [7] has two disadvantages. First, this equation is used in connection with 
the function defined in Appendix B of [7]. However, is only well defined if i/o is not 
a root of the argument of the second logarithm. In fact, there are several cases where the 
second logarithm of S3 is divergent but the final result of T is finite. This point is made 
clear from the calculation of the function Si in Appendix [Bj Second, this equation leads 
to the evaluation of several 77-functions which might suffer from nasty cancellations. In 
particular, the 77-functions in eqs. (6.26, B.2) of [7J may cause serious problems. 

The representation ([H]) leads naturally to two ways of numerical evaluation. This will 
be discussed in Section |3l 

2.3.2 Special cases: A = or B = 

These special cases can occur when we apply Eq. ( 1T0|) for the case of more than one 
lightlike momenta. 

If i? = then both denominators in Eq. (!24l) are linear in y then the ?/-integral can 
be easily performed. The result reads 

nAO,C;G,H;D,E,F,J) = _ ^g^^.^^ _ E ^-l^ 



1 

X / dy ln(A,y2 + 5,y + C,), (30) 

^0 y ~ Hi 



where are the two roots of 

{HA - CG)y^ + {HD -CJ - EG)y + H{F - ie) - E{J - le) = (31) 

or 

-{HD -CJ- EG) - ^/det{Q) _ -{H D - C J - EG) + 
~ 2{HA - CG) ' ~ 2{HA - GG) ' 

det(g) = {HD -GJ - EGf -A{HA-GG){HF - EJ). (32) 

The discriminant is the determinant of the Landau matrix Q if our calculation is set in 
the context of Eq. (llOp or Eq. (1191) . The other coefficients are given in Table |H 

If y4 = 0, we change the integration order by using Eq. f|T3|) and then change the 
integration variables in two successive steps: x = x' + y and y = 1 — y'. At the end, one 
should obtain the following relation 

r(0, B, G, G, H;D, E, F, J) = r{B + C, 0, -C, -G - H, G; 

- D - E -2B -2G,D + G,D + E + B + G + F,G + H + J). ^ ' 

Thus, we get back to the above case of i? = 0. 
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Table 4: The coefficients of the T function 



j 


A, 


5, 


G, 


1 


A + C 


E + D 


F-ie 


2 





G + H 


J -ie' 


3 





G 


J-ie' 


4 


A 


D 


F-ie 



We conclude this section by observing that the above results can be also applied for 
the case of real masses. In this case, the two poles yi^2 in the rhs of eqs. f|T^ [25| 1301) can 
be real and the ze-prescription is needed. The position of the poles depend on the value of 
internal masses and external momenta. Thus, it might happen that the two poles become 
equal and pinch the integration contour, [0, 1]. If this occurs then the Dq function can 
have a four-point leading Landau singularity (LLS). The issue of Landau singularities will 
be discussed in the appendices. 

3 The code DOC 

The above method has been implemented into a Fortran 77 code, DOC, which can calcu- 
late the scalar one-loop four-point integral, Dq, with complex/real internal masses. The 
restriction is that at least one external momentum must be lightlike or a must be real (see 
Eq. (l2Ti) ). The code is best exploited when included into the library LoopTools where a 
reduction procedure to calculate tensor one-loop five-point integrals with complex masses 
is available. Thus, the new version of LoopTools (LoopTools-2 . 4) can handle all the 
one- loop integrals with complex masses up to five-point. The original code DOC can be 
downloaded at: 

|http://wwwth. mppmu.mpg.de/members/ldninh/index.html 
LoopTools-2.4 including DOcl§ can be downloaded at: 

http:// www, feynarts. de/looptools/l 

In the language of LoopTools, the Dq function with complex masses, also named DOC, 
is invoked in the same way as the usual DO with real masses. However, all the arguments 
of DOC must be of double complex type. 

The implementation of Eq. (l25l) is straightforward. However, there is an important 
point we want to stress. The residue term on the rhs may cause serious problem of 
numerical cancellation when the imaginary part of the i?jS become very small. Thus we 
want to eliminate this term as far as possible. Our algorithm is based on the simple fact 
that this contribution vanishes if the pole is outside of the beta triangle. If (3 is real, the 
triangle becomes a line in the complex plane. The residue contribution vanishes. If j3 is 
complex, there are two choices. First, we take j3 = jSi which is the largest root of Eq. fl28l) . 

^Thomas Hahn has adapted this code to LoopTools. The structure is unchanged. The way of coding 
is shghtly modified. The numerical results of the two codes should be identical. 
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There is a very efficient method to check whether a point (pole) is inside a triangle by 
using the barycentric technique. This idea is very simple. From a triangle, ABC, one can 
define a barycentric reference frame {AB, AC) where A is the origin. The coordinates of 
a point P, {x, y), are easily calculated by writing: 

AP = xAB + yAC. (34) 

Multiply both sides by AB and AC to get a system of two equations to solve for (x, y). 
The conditions for P to be inside the triangle ABC are: x and y must be positive and 
X + y < 1. However, one sees that the problem of numerical cancellation is still present 
if x or y is very small, i.e. the point P is very close to the border of the triangle. If this 
happens, we repeat the above steps with the second choice P = ^2- That is what we have 
done in the code. 

In the worst case, if x or y is very small for both choices of f3 then the code will produce 
a warning of numerical cancellation but still continue the computation. The result in this 
case is not reliable and must be cross-checked. The code includes a means to do this. 
Indeed, we have written two different codes, called versions 1 and 2, to evaluate the rhs 
of Eq. (1251) based on two different strategies. The first one is described as above. In the 
second one, the residue contribution is always calculated as long as (3 is complex. One 
has to specify which version when calling the main function: 

d0c=D0Cn(meth,pl,p2,p3,p4,plp2,p2p3,ml,m2,m3,m4) 

where meth=0 (version 1) or 1 (version 2). The two versions are different only when (3 is 
complex. For LoopTools-2 . 4, one can choose different versions by defining a version key 
before calling the DOC function as 

call setversionkey(i*KeyDOC) 
d0c=D0C(pl,p2,p3,p4,plp2,p2p3,ml,m2,m3,m4) 

where i can be (compute version 1), 1 (compute version 2), 2 (compute both versions, 
return 1) or 3 (compute both versions, return 2). The default value is i=0. 

The structure of the original code, written in double precision Fortran 77, is as follows. 

• DOCn.F: the main function to return the value of the integral Dq as a function of 
masses and momenta (see Eq. ([1])). It rotates the momenta to the right position and 
prepares input parameters for Eq. ([5]). 

• dOc_abc.F: a subroutine called in DOCn.F to calculate the rhs of Eq. 

• fT13.F is a subroutine to calculate the function D^^^ as defined by Eq. (11 ip . 

• fTn.F (including fTn_lin.F to handle special cases) is a subroutine to calculate the 
T function (defined by Eq. (l23l) ) by using the barycentric technique described above. 
It belongs to version 1. 

• fTn_resd.F(including fTn_lin.F to handle special cases) is another subroutine to 
calculate the T function for cross-checking purpose. It belongs to version 2. 

• fSSn.F and fR.F are subroutines to calculate the functions Si and TZi defined in 
Appendix [Bl 
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• f S2 . F is a subroutine to calculate the functions T and Dq if the integrand is just 
a logarithm, i.e. the prefactor is just a constant. 

• Fundamental functions including Spence function, logarithm, ?7-function are taken 
from FF. 

The code has been carefully checked as follows. Since the ze-prescription is kept explicitly 
in our code, the code works also for the case of real masses. In this case, our results have 
been checked by comparing with the DO function of LoopTools. Full agreement with the 
version 'b' of LoopTools has been founcj^. We then introduced the widths for internal 
particles and checked the results in two limits: very small and very large widths. In the 
former case, we recovered the results of real masses. For the latter case, the results were 
checked against the ones of a routine where the integration is done numerically by using 
BASES [13]. We observed good agreement within the integration error (see Tabled]). 

Tables l^fTfSl give some numerical results of the code for two cases: complex and real 
masses. For the latter we observed that our code agrees with version 'b' of LoopTools 
while disagrees with version 'a'. In fact, our code produces the same results with both 
versions of LoopTools in almost all physical cases. The results shown in Tables ITfHl are 
a few cases where disagreement has been founcf]. Thus, care must be taken when using 
loop libraries and one should always cross-check the results by using different methods. 
For instance, it can be very useful to check the results of real masses in a tricky situation 
by using the code for complex masses with a very small imaginary part. 

Table 5: Comparisons with a numerical integration routine using BASES, complex masses 
with large imaginary parts: mj = 28900-25500i, = 6400-16000i, m| = 8100-22500^, 
ml = 40000 - SOOOOi, pj = 100, pi = 32400, pi = 400 pj = 250000, {p2 + p^f = 84100. 
Numerical integration errors are about (0.03%, 0.07%). 



(Pl +P2f 


DOC 


Numerical integration 


48400 


(-1.50823957E-10,1.62902558E-11) 


(-1.5077496E-10,1.63063608E-11) 


52900 


(-1.51384299E-10,1.37108163E-11) 


(-1.51339092E-10,1.37149562E-11) 


57600 


(-1.51834347E-10,1.09709949E-11) 


(-1.51787895E-10,1.09683749E-11) 


62500 


(-1.52150627E-10,8.07814957E-12 ) 


(-1.52101762E-10,8.07660917E-12) 



^LoopTools has two versions, 'a' and 'b', to calculate the DO function [TT]. Version 'a' is based on FF 
while version 'b' uses the method of [H]. 

"LoopTools does produce some warning or error messages in these cases. 
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Table 6: Some results of DOC, complex masses with small imaginary parts: ml = 28900 — 
255i, ml = 6400 - 160i, mj = 8100 - 225i, mj = 40000 - 800z, pj = 250000, {pi +p2f = 
72900, ip2+P3y = 84100 



p'i 


pI 


p'i 


DOCxlO-^° 


100 


32400 


400 


(-37.6230835,-134.501524) 





32400 


400 


(-37.8749852,-134.892905) 





32400 





(-38.6796394,-136.041222) 











(-28.9137581,16.3570695) 



Table 7: Comparisons between LoopTools and DOC, real masses: mf = m\ = m\ = ml = 
30276, pI = 14400, (^1+^2)^ = 57600 





(0,124609,20,129630) 


(0,0,20,73473) 


LoopTools 'a'xlO"^" 


(-8.4106454,227.290811) 


(19.8584112,-8.23191831) 


LoopTools 'b'xlO^^'^ 


(-8.4106454,26.7103536) 


(3.32391023,-6.755152E-16) 



000x10'^*^ (-8.4106454,26.7103536) (3.32391023,-2.067952E-15)^ 



Table 8: Continuation of Table [71 



{pi,pi,pi,{p2 +P3)'^) 


(0,0,0,734730) 


(100,50,20,734730) 


LoopTools 'a'xlO^^'J 


(-0.400586866,-6.57420392) 


(-0.401370928,-6.57507623) 


LoopTools 'b'xlO"^" 


(-0.400586866,1.85027075) 


(-0.401370928,1.85083346) 


000x10"^*^ 


(-0.400586866,1.85027075) 


(-0.401370928,1.85083346) 



4 Conclusions 

The Fortran code DOC allows to calculate the scalar one-loop four-point function with 
complex masses for all physical processes at colliders. The method and structure of the 
code were explained. The results related to massless external particles and the discussion 
of Landau singularities complement the study of 't Hooft and Veltman. The code has 
been included in LoopTools-2 .4. Please refer to this paper if you use the code. Bug 
reports or comments should be sent to the author^. 
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Appendices 



A Four-point leading Landau singularity 

The method of this paper can also be applied for the case of real internal masses. In 
this case the ie-prescription is needed to define various multivalued functions (logarithms, 
Spence functions) in the final results. From the practical viewpoint, this method is not 
as efficient as other methods based on projective transformation [3 |T2] since the results 
are written in terms of a larger number of Spence functions. However, it has some use- 
ful virtues as a direct calculation and we will show that this method is indeed a very 
easy means to understand the analytical properties of the scalar four-point function. To 
illustrate this purpose, we first consider the following function 



^0 Ky"^ + Ly + N -ie' 
An important property of this function is that the denominator has, in general, two poles 



-L - VL^ - 4:KN 
2/1 = 



2K ReVU^lKN 
-L + ^/iX^WN le' 
= ^ ^ ReVL'^-AKN ^^''^ 

If these poles are real then it may occur that, in some region of parameter space, they pinch 
the integration contour JO, 1], leading to a singularity. This phenomenon occurs when the 
two roots become equaH, i.e. det(Q) = — AKN — > 0. Since y2 — Vi oc 2ie' ^ 0, we can 
split T\ into two terms as usual 

Ti = — -i ^ C dy (— ln( V + By + C-ie). (A.3) 

^(2/2-2/1) Jo V2/-2/2 2/ -2/1/ 

Using the results of Appendix [Bl and taking the limit 2/2^2/1 = 2/0 = —L/{2K), yo is 
real, we get 



rr'- = ^rr^'' MA - le) + r/(-^i, -z,) + M{z,, y,) + M{z,, y,)] (A.4) 
^det(Q) 



where ei = e'/ReV-^^ — AKN, 2:1 2 are the roots of the argument of the logarithm and 

[ - ln(-«e^) - 7r«sign(e^) it z = yo 
r/(a,6,c) = 27ri[^(-Ima)^(-Im6)^(Imc) -^(Ima)^(Im6)^(-Imc)]. (A.5) 



WThis cannot happen in the case of complex masses since the discriminant, the Landau determinant 
det((5), is protected from zero by an imaginary part. 
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We remark that if 

yoivo - 1) < 0, or < Vo < 1 (A.6) 

then the result is divergent. This is the origin of the four-point LLS. We also observe that 
if yo = z\ or/and then we have the coincident of four-point and three-point Landau 
singularities. The latter will be discussed in Appendix [Bl 

The scalar one- loop four-point function, Dq? is a sum of T^s. The above conditions 
for leading Landau singularity are therefore necessary but not sufficient. In fact, the 
conditions for scalar one-loop integrals to have a LLS can be obtained by using Landau 
equations [H [151 E] • However, one might expect to find out the same conditions from a 
direct calculation. To illustrate the key points while keeping our explanation as simple as 
possible, we consider the case of Dq function with two opposite lightlike external momenta, 
see Eq. (IT^ . We ignore the case of coincident singularity and assume y^ ^ Z\^2-, Eq- (IA.4P 
becomes 



rr'- = rr^ ' HAy'o + Byo + C- ze), (A.7) 

^/det{Q) 

where we recall that yo is real. Apply this result to Eq. f[T^ we get 



° ^/d^t{Q) V Uo + Vo-te Uoyo + Vo-ieJ 

where the index in the rhs means that the functions are real (without ie) and calculated 
at y = yo. We now use the following property 

to get 

° V^d^ 
Thus, necessary and sufficient conditions for the function -Dg^^'' to have a LLS are 

det(Q) = 

< ?/o < 1 (A.ll) 
(f/o + K))(t/ol/o + K)) <0 and ToK, < 

where yo is the root of the equation SV — TU = 0. We have checked that this result 
agrees with the general conditions obtained by using Landau equations |6j. 



B Basic integral and three-point Landau singularity 



We calculate the following integral 

y-yo + ie' 
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where A is real, while B, C and may be complex, with the restriction that the imaginary 
part of the argument of the logarithm has always the same sign for < ?/ < 1, as indicated 
by the " — ie" . The " — ie'" is only needed when uq becomes real. 
The numerator can be written as 

\n{Ay'^ + By + C -ie) = \n{A - ie) + ln{y - zi) + \n{y - Z2) + r?(--2i, -Z2), (B.13) 

where zi^2 are the two roots of the argument of the logarithm. We have 

Si = [\n{A-ie)+r]{-zi,-Z2)]\n^^^^ + ni{yo,zi)+ni{yo,Z2), (B.14) 
where the T^i-function is defined as 

TZi{y,z)= dx^ — ^ B.15 

Jo x-y + iey 

with ez^y are infinitesimal, li y = z and being complex, we use 

Uiiy, z) = ^ln (^^^ [ln(l - z) + In(-z)]. (B.16) 
li y = z and being real, we then use 

ln(a; — z + iez) = ln(a; — z + iey) — r]{x — z + iez, : — ) (B-17) 



to get 



I / 1 — y \ ie 



1€ 



:b.i8) 



We observe that iiO<y = z<l and tyCz < then integration contour is pinched and 
this singularity corresponds to the anomalous threshold of scalar three-point integral (also 
called the three-point Landau singularity) whose nature is logarithmic. Since the value 
of the 77-function is purely imaginary, the imaginary part of TZi is divergent while there is 
a finite jump in the real part of TZi at the singular point. 
li z y then the result reads 

TZi{y, z) = Hy - z) In - + 7^(l/, z) (B.19) 

where the function TZ{y,z) is given in Appendix B of [7j. When calculating the rhs of 
Eq. flB.191) . if the argument {y — z) becomes negatively real then one has to introduce an 
infinitesimal imaginary part ieyz- The result does not depend on the sign of this imaginary 
part as long as we introduce it everywhere in the rhs of Eq. flB.191) . Thus, the simplest 
way is to write y — z = y — z + iez- This completes our calculation of Si. 
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